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ABSTRACT 

We investigate the rms peculiar velocity of galaxy clusters in the Lambda cold dark 
matter (ACDM) and tau cold dark matter (tCDM) cosmological models using N-body 
simulations. Cluster velocities for different cluster masses and radii are examined. 
To identify clusters in the simulations we use two methods: the standard friends-of- 
friends (FOF) method and the method, where the clusters are defined as the maxima 
of the density field smoothed on the scale R ~ lh^ 1 Mpc (DENSMAX). If we use 
the DENSMAX method, the size of the selected clusters is similar for all clusters. 
We find that the rms velocity of clusters defined with the DENSMAX method is 
almost independent of the cluster density and similar to the linear theory expectations. 
The rms velocity of FOF clusters decreases with the cluster mass and radius. In 
the ACDM model, the rms peculiar velocity of massive clusters with an intercluster 
separation d c i = 50/i _1 Mpc is «15% smaller than the rms velocity of the clusters 
with a separation d c i = 10/i -1 Mpc. 
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1 INTRODUCTION 

One of the interesting unknowns in cosmology is the large- 
scale peculiar velocity field in the Universe. The peculiar 
velocity field can be studied by using galaxies or clusters of 
galaxies. However, there are some advantages in studying 
the peculiar velocity field by using galaxy clusters. One of 
these advantages comes from the fact that, on scales probed 
by galaxy clusters, velocity fluctuations are largely in the 
quasi-linear regime and close to the initial state from which 
large scale structures developed. In addition, peculiar ve- 
locities of clusters can be determined more accurately than 
peculiar velocities of galaxies since the distance to each clus- 
ter can be obtained from a large number of member galaxies, 
thus considerably reducing the velocity uncertainties of clus- 
ters. Cluster motions could therefore provide an important 
tool in probing the large-scale peculiar velocity field. 

Peculiar velocities of clusters of galaxies have been stud- 
ied in several papers (e.g. Bahcall, Gramann & Cen 1994; 
Lauer & Postman 1994; Bahcall & Oh 1996; Moscardini et 
al. 1996; Borgani et al. 1997; Watkins 1997; Dale et al. 1999; 
Hudson et al. 1999; Borgani et al. 2000; Colless et al. 2001). 
Watkins (1997) developed a likelihood method for estimat- 
ing the rms peculiar velocity of clusters from line-of sight 
velocity measurements. This method was applied to two ob- 
served samples of cluster peculiar velocities: the SCI sample 
(Giovanelli et al. 1997) and a subsample of the Mark III 



catalogue (Willick et al. 1997). Watkins (1997) found that 
the rms one-dimensional cluster peculiar velocity is 256^^ 
km s , which corresponds to the three-dimensional rms ve- 
locity 459tij^ km s_1 ■ Dale et aL (1999) obtained Tully- 
Fisher peculiar velocities for 52 Abell clusters distributed 
over the whole the sky between ~ 50 and ~ 200/i _1 Mpc. 
They found that the rms one-dimensional cluster peculiar 
velocity is 341 ± 93 km s^ 1 , which corresponds to the three- 
dimensional rms velocity 591 ± 161 km s _1 . 

Radial peculiar velocities of clusters can be determined 
to large distances by measuring the Sunyaev-Zeldovich 
(1980) (SZ) effect. Rephaeli & Lahav (1991) made one of 
the first estimates of the possibility of measuring the pecu- 
liar velocities by using the SZ effect for a selected sample of 
galaxy clusters. However, most convincing measurements for 
individual clusters have been done only recently, using the 
new generation of sensitive bolometers (Holzapfel et al. 1997; 
Lamarre et al. 1998). The accuracy of SZ measurements for 
determining peculiar velocities of clusters have been stud- 
ied in several papers (Haehnelt & Tegmark 1996; Aghanim, 
Gorski & Puget 2001; Diego et al. 2002; Holder 2002; Na- 
gai, Kravtsov & Kosowsky 2002). With fiK sensitivity on 
arcminute scales at several frequencies it will be possible to 
measure peculiar velocities to an accuracy of ~130 km s . 

In this paper we study the rms peculiar velocity of clus- 
ters, v rm3 , in different cosmological models assuming that 
the initial density fluctuation field is a Gaussian field. To in- 
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vestigate the nonlinear regime, we use N-body simulations. 
We examine cluster peculiar velocities for different cluster 
masses. Do cluster velocities depend on their masses? The 
rms peculiar velocity of peaks in the initial Gaussian field 
does not depend on the height of peaks (Bardeen et al. 1986). 
However, as an initially Gaussian density field evolves gravi- 
tationally in the nonlinear regime, it becomes non-Gaussian. 

The evolution of peculiar velocities of galaxy clusters 
in different N-body models with an initially Gaussian den- 
sity field has been examined in several papers (e.g. Bahcall, 
Gramann, Cen 1994; Croft & Efstathiou 1994; Suhhonenko 
& Gramann 1999, Colberg et al 2000; Sheth & Diaferio 
2001). Groth & Efstathiou (1994) studied the cluster pecu- 
liar velocities for two different cluster richnesses, described 
by the mean intercluster separations d c i = 30fe _1 Mpc and 
d c i = 55/1" 1 Mpc. They found that the peculiar velocity dis- 
tributions are almost independent of cluster richness. Sheth 
& Diaferio (2001) studied rms peculiar velocity of clusters 
for different masses and found that in N-body simulations 
the rms velocity of clusters depends weakly on cluster mass, 
with a small trend that for more massive clusters the rms 
peculiar velocity decreases. On the other hand, Suhhonenko 
& Gramann (1999) investigated the properties of clusters 
using N-body simulations and found that the rms peculiar 
velocities of clusters increase with cluster richness. 

In this paper we study the dependence of v rm3 on clus- 
ter masses in more detail. We also examine the rms peculiar 
velocities of clusters for different cluster radii. We use the 
N-body simulations published by the Virgo Consortium and 
described in detail by Jenkins et al. (1998). These simu- 
lations were carried out using a parallel, adaptive particle- 
particle/particle-mesh (AP 3 M) code (Couchman, Thomas & 
Pearce 1995; Pearce & Couchman 1997). In this paper we 
analyze the velocities in the ACDM and rCDM model (see 
Jenkins et al. (1998) and Section 2 for the description of 
the cosmological parameters in these models). We also fol- 
low the evolution of particles in a similar ACDM model with 
the same cosmological parameters but using a particle-mesh 
(PM) code. In this PM code we use a traditional two-point 
finite-difference approximation to calculate the forces on the 
grid. We also examine the shifted-mesh scheme (see eq. (13) 
below) and show that this scheme artificially boosts the ve- 
locities of clusters. 

To identify clusters in the simulations we use two meth- 
ods: the standard friends-of-friends (FOF) method and the 
method, were the clusters are defined as the maxima of the 
density field smoothed on the scale R ~ lh~ l Mpc (DENS- 
MAX). To determine the velocity of DENSMAX clusters, wc 
use the same smoothing scale as for the density field. The 
velocity of FOF clusters is defined to be the mean velocity 
of all the particles in the cluster. 

This paper is organized as follows. In Section 2 we 
describe the cosmological models, N-body simulations and 
cluster selection algorithms used. We also examine the linear 
theory predictions for the peculiar velocities of peaks in the 
Gaussian field. In Section 3 we examine the rms peculiar 
velocity of clusters for different cluster masses and radii, 
and compare the cluster velocities with the linear theory 
predictions. In Section 4 we briefly discuss the effect of the 
shifted-mesh scheme on the cluster velocities. A summary 
and discussion are presented in Section 5. 



2 MODELS 
2.1 Simulations 

We analyze peculiar velocities in N-body simulations carried 
out by the Virgo consortium for two cosmological models, 
ACDM (fi = 0.3, A = 0.7) and rCDM (fi = 1) as de- 
scribed by Jenkins et al. (1998). In these cold dark matter 
(CDM) models the power spectrum of the initial conditions 
was chosen to be in the form given by Bond & Efstathiou 
(1984), 



V ' [l + iaq+ibq^^ + icq) 2 )-] 2 /^ V ; 

where q = k/T, a = 6.4ft" 1 Mpc, b = 3/T 1 Mpc, c = 1.7ft- 1 
Mpc, v = 1.13 and Y = 0.21. The normalization constant, 
A, was chosen by fixing the value of as (the linearly ex- 
trapolated mass fluctuation in spheres of radius 8ft _1 Mpc); 
a s = 0.9 and a s = 0.51 for the ACDM and rCDM model, 
respectively. 

We investigated the linear theory predictions for pecu- 
liar velocities of peaks in the ACDM and rCDM model. The 
linear rms velocity fluctuation on a given scale R can be ex- 
pressed as 

a v (R) =flb/(no)cT-i(Ji), (2) 

where the spectral moments aj are defined for any integer j 
by 

a) = ^_ J P(k)W 2 (kR)k 2j+2 dk, (3) 

W(kR) is a window function and /(fio) is the dimensionless 
growth rate. The function /(fio) = 0.51 and /(fio) = 1.0 in 
the ACDM and rCDM model, respectively. (We note that 
the approximation /(fio) = fio' 6 underestimates the dimen- 
sionless growth rate by ~ 5% in the flat fi = 0.3 model). 

Bardeen et al. (1986) showed that in the linear approx- 
imation the rms peculiar velocity at peaks of the smoothed 
density field differs systematically from a v (R), and can be 
expressed as 

a p (R) = a v (R) ^\-atla\a^. (4) 

In this approximation, the rms velocities of peaks do not 
depend on the height of the peaks. 

Fig. 1 shows the rms peculiar velocity of peaks, a p (R), 
for the ACDM and rCDM model. We have used the top- 
hat window function. For comparison, we show also the rms 
peculiar velocity a v (R) for the same models. For the radius 
R — lft -1 Mpc, a p = 509 km/s and a p — 562 km/s in the 
ACDM and rCDM models, respectively. At this radius, a v 
is lower than a v about ~ 2 per cent for the models studied. 
On larger scales, the difference between a p and a v increases. 

Next we study peculiar velocities in N-body simulations. 
The Virgo simulations were created using a parallel adaptive 
particle-particle/particle-mesh (AP 3 M) code as described by 
Couchman, Thomas & Pearce (1995) and Pearce & Couch- 
man (1997). It supplements the standard P M algorithm 
(Efstathiou et al. 1985) by recursively placing higher resolu- 
tion meshes, 'refinements', in heavily clustered regions. The 
Virgo simulations were done on two large Cray T3D par- 
allel supercomputers at the Edinburgh Parallel Computing 
Center and at the Computing Center of the Max Plank Soci- 
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Figure 1. The rms peculiar velocity of peaks, a p (R), in the 
ACDM model with Co = 0.3 and as = 0.9 (heavy solid line) 
and in the tCDM model with Qo = 1.0 and erg = 0.51 (heavy 
dot-dashed line). The light curves show the corresponding rms 
peculiar velocity cr v (R) for the same models. 



ety in Garching. These simulations are publicly available at 



http :/ /www. mpa- garching. mpg. de/ Virgo/ virgoproject. htm . 

In the simulations used here, the evolution of particles 
was followed in the comoving box of size L = 239. 5/i -1 Mpc. 
The number of particles was N p = 256 3 . Therefore, the mean 
particle separation X p — L/Np^ a = 0.9355/i _1 Mpc. The 
mass of the particle m p — pb\p = 6.82 x 10 10 ft -1 Mq and 
m p = 2.27 x ICF/i^Mq in the ACDM and rCDM mod- 
els, respectively (here pt is the mean background density). 
The gravitational softening length is r so f t = 25h~ 1 kpc and 
r soft = 36/i _1 kpc , respectively (see Jenkins et al. (1998) for 
a detailed description of the force calculation scheme used 
in the Virgo simulations). We denote the Virgo ACDM and 
rCDM models as the model ACDM1 and rCDM. 

We also investigated the evolution of 256 3 particles on 
a 256 3 grid using a particle-mesh (PM) code described by 
Gramann (1988) and Suhhonenko & Gramann (1999). The 
PM code achieves the force resolution close to the mean 
particle separation X p . The PM method is discussed in detail 
by Hockney & Eastwood (1981) and Efstathiou et al.(1985). 
The cosmological parameters for this PM simulation were 
chosen similar to the ACDM1 model. We chose the flat fio = 
0.3 model and used the initial power spectrum P(k) given 
in eq. (1) (with ug = 0.9). The comoving box size was L = 
239. 5/i." 1 Mpc. Therefore, the mean particle separation, \ p , 
and the mass of the particle, m p , in this model are the same 
as used in the ACDM1 model. We denote this model as the 
ACDM2 model. 

We examined the rms velocity of particles in the sim- 
ulations studied. In the ACDM1 and rCDM models, the 
rms velocity of particles was 648 kms -1 and 636 kms -1 , 
respectively. In the ACDM2 model, the rms velocity was 
575 kms -1 . Due to the small-scale smoothing inherent to 
the PM method, the intrinsic velocity dispersions of clusters 
in the ACDM2 model are smaller than in the ACDM1 model. 
The velocity field of clusters in these models is expected to 
be similar. 

In the ACDM2 model we used the traditional two-point 
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Figure 2. The cluster mass function in the Virgo ACDM model 
for the parameter b = 0.164 (solid line) and in the rCDM model 
for b = 0.2 (dot-dashed line). For comparison we show the fitting 
formulae given by Jenkins et al. (2001) for the mass function in 
the ACDM (dashed line) and in the tCDM model (dotted line). 



approximation to calculate the acceleration on the grid. We 
also studied the evolution of particles using a PM code 
with the shifted-mesh scheme. The effect of the shifted-mesh 
scheme on the cluster velocities is discussed in Section 4. 



2.2 Selection of clusters 

We used two different algorithms to identify clusters in sim- 
ulations: the standard friends-of-friends (FOF) algorithm, 
and the algorithm, where clusters are defined as maxima 
of the density field smoothed on the scale R ~ 1/tA 1 Mpc 
(DENSMAX). 

The friends-of-friends group finder algorithm was ap- 
plied using the program suite developed by the cosmology 
group in th e University of Washington. These prog rams are 
available at 



http : / '/www-hpec. astro. Washington, edv 



The FOF cluster finder depends on one parameter b, 
which defines the linking length as b\ p . The conventional 
choice for this parameter is b = 0.2 (see e.g. Gotz, Huchra & 
Brandenberger 1998; Jenkins et al. 2001). In this paper we 
also define clusters by using the value b = 0.2. We also study 
velocities of the clusters defined by the parameters b — 0.15 
and b = 0.3. In the limit of very large numbers of particles 
per object, FOF approximately selects the matter enclosed 
by an isodensity contour at 1/fe 3 . 

To test our FOF output data, we found the mass func- 
tion of clusters in the Virgo simulations. The cluster mass 
function in these simulations has been studied in detail by 
Jenkins et al. (2001). Fig. 2 shows the mass function of 
clusters determined for b — 0.164 in the ACDM model and 
for b = 0.2 in the rCDM model. For comparison we show 
the mass functions given by the approximations obtained 
by Jenkins et al. (2001) (eq. (B2) and (Bl) for the ACDM 
and rCDM model, respectively). Jenkins et al. (2001) stud- 
ied the mass function at the high mass end up to the point 
where the predicted Poisson abundance errors reach 10%. 
In the simulations studied here, this limit is reached when 
the number density of clusters is n = 7.3 x 10" 6 ft 3 Mpc- 3 . 
Fig. 2 shows that if n is larger than this value, the agree- 
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ment between our results and these obtained by Jenkins et 
al. (2001) is very good. 

We studied clusters that contained at least ten particles. 
The three-dimensional velocity of each cluster was defined 

as 

i=i 

where N p is the number of particles in the cluster and Vi is 
the velocity of the particle i in the cluster. To characterize 
the size of the cluster, we use the effective radius defined as 

R *ff = - ^ 2 + (»* - + (* - *) 2 ] 1/2 > ^ 

p i=l 

where Xi,yi, Zi are the particle coordinates in the cluster and 
(x, y, z) are the coordinates of the cluster centre. If we use 
the FOF method, the mean size for the high-mass clusters 
is larger than the mean size for the low-mass clusters. 

We also selected clusters using the DENSMAX method. 
In this case clusters were identified in the simulations as 
maxima of the density field that was determined on a 256 3 
grid using the cloud-in-cell (CIC) scheme. To determine pe- 
culiar velocities of clusters, v c i, we calculated the peculiar 
velocity field on a 256 3 grid using the CIC-scheme, and found 
peculiar velocities at the grid points where the clusters had 
been identified. If we use the DENSMAX method, the size 
of the selected clusters is similar for all clusters, given by 
the cell size, X g . In our case, the cell size X g = 0.9355/i _1 
Mpc. 

To determine the rms peculiar velocities of clusters, we 
used the equation 

2 2,2 1 V"^ 2,2 , n \ 

V rm s=V s +V L = —^V cli +V L , (7) 

i= 1 

where the parameter v s describes the dispersion of cluster 
velocities, v c u, derived from simulation and the parameter 
dl is the linear contribution from velocity fluctuations on 
scales greater than the size of the simulation box L and is 
given by 

N c i is the number of clusters studied. By using eq. (4), the 
linear rms peculiar velocity of peaks can be written as 

a 2 p (R)=al(R)-Hlf 2 (n )^§-. (9) 

The second term in this expression is not sensitive to the 
amplitude of large-scale fluctuations at wavenumbers k < 
2-k/L. Therefore, the linear rms velocity of peaks can be 
expressed as 

al(R) = a' p 2 (R) + vl, (10) 

where <j' p (R) is determined by the power spectrum at 
wavenumbers k > 2n/L and vl is given by eq. (8). For 
ACDM and rCDM models, we found that v L = 220 kms -1 
and vl — 245 kms -1 , respectively. 

If the one- dimensional velocities of clusters, v X i, follow 



Table 1. The number of clusters, N c i, in different density and 
mass intervals. iV c ; for different density intervals is given for the 
clusters selected with the DENSMAX method. N cl for different 
mass intervals is given for the clusters determined with the FOF 
method with b = 0.2. 
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10 14 - 5 x 10 14 
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5 x 10 14 - 10 15 
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18 


24 


10 15 - 5 x 10 15 
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5 


3 



a Gaussian distribution with a mean v x = and a dispersion 
a 2 , then the sum 

* = v2 " 

i=l 

is distributed as a \ 2 distribution with the number of degrees 
of freedom v = 3iV c (. In this case, the rms error for the 
variable v 2 ms can be determined as 

^-^{m:^ 1 (12) 

As a first step, the one-dimensional distribution of cluster 
velocities can be approximated as a Gaussian distribution 
(see e.g. Bahcall, Gramann, Cen (1994) for the study of 
the velocity distribution of clusters in different cosmologi- 
cal models). We used eq. (12) to estimate the error bars for 
the rms velocities of clusters. 



3 RESULTS 

First, we investigated the rms velocity of clusters in different 
density and mass intervals. The results are presented in Fig. 
3. The rms velocities are shown for the intervals, where the 
number of clusters is N c i > 10. 

The clusters defined with the DENSMAX method were 
divided into subgroups according to their density. We stud- 
ied the rms velocity of clusters in seven subgroups, where the 
density p/pb was in the range 1 — 5, 5—10, 10 — 50,...,1000 — 
5000. Table 1 shows the number of clusters and the upper 
right panel in Fig. 3 shows the rms velocity of clusters in 
different density intervals. We see that the rms velocities of 
clusters in the ACDM1 and ACDM2 models are similar. In 
the models studied, the rms velocity of clusters is almost 
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Figure 3. The rms peculiar velocities of clusters for different densities and masses in the ACDM1 model (solid lines), tCDM model 
(dot-dashed lines) and in the ACDM2 model (dashed lines). The upper left panel shows the results for the clusters defined by the 
DENSMAX method and the lower left panel shows the results for the clusters selected by the FOF method with b = 0.2. The right 
panels show the rms velocities of clusters defined by the FOF method with b = 0.15 and b = 0.3, respectively. 



independent of the density. The rms velocity somewhat in- 
creases at smaller densities. However, this increase is small 
(fs 10 per cent). For the range p/pb = 100 — 500, the rms 
velocity is 505 km s _1 and 570 km s _1 in the ACDM1 and 
rCDM models, respectively. These values are similar to the 
linear theory expectations. At the radius R — lh^ 1 Mpc, 
the rms peculiar velocity of peaks a p = 509 km s _1 and 562 
km s _1 in the ACDM1 and rCDM models, respectively. The 
rms velocity for the low-density clusters is somewhat smaller 
than predicted by the linear theory. 

The clusters determined with the FOF method were 
divided into subgroups according to their mass. We studied 
the rms velocity of clusters in eight subgroups, where the 
mass was in the range (5 x 10 11 - 10 12 )/i _1 M o , (10 15 - 
5 x lO 15 )ft~ 1 M0. Table 3 shows the number of clusters and 
Fig. 3 demonstrates the rms peculiar velocity of clusters in 
different mass intervals. The lower left panel in Fig. 3 shows 
the results for the clusters determined by b — 0.2. The right 
panels show the rms velocities of clusters defined by b = 0.15 
and b = 0.3. 

We see that in the ACDM1 and ACDM2 models the 
rms velocities are similar. The rms velocity of FOF clusters 
decreases with cluster mass. The rms velocity of massive 



FOF clusters is smaller than the rms velocity of high-density 
clusters determined with the DENSMAX method. In the 
ACDM model for b — 0.2, the rms velocity of clusters is 
525 kms -1 in the mass interval (5 x 10 11 - W 12 )h' 1 M Q and 
430 kms" 1 in the mass interval (10 15 - 5 x 10 15 )/i _1 M o . 
For b — 0.15 and b = 0.3, this effect is similar. These results 
are in good agreement with the results obtained by Sheth & 
Diaferio (2001). They studied the rms velocities of clusters 
in different mass intervals and found that the rms cluster 
velocity decreases with mass. 

We also studied the rms velocity of clusters with den- 
sity (or mass) higher than a given threshold density (or 
mass). The results are presented in Fig. 4. The DENS- 
MAX clusters were ranked according to their density, and 
we selected N c i = (L/d c i) highest ranked clusters to pro- 
duce cluster catalogues with a mean intercluster separation 
10 — 50fe _1 Mpc. Similarly, the FOF clusters were ranked ac- 
cording to their mass. Table 2 shows the density and mass 
thresholds used to produce cluster catalogues for different 
values of the mean cluster separation. For comparison, the 
number density of observed APM clusters and Abell clus- 
ters is n ct ~ 3.4 x 10~ 5 /i 3 Mpc~ 3 (d ct ~ 31/i _1 Mpc) and 
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Figure 4. The rms peculiar velocities of clusters for different values of the mean cluster separation, d c i . The clusters arc ranked according 
to their density (DENSMAX clusters) or mass (FOF clusters). The lines are defined as in Fig. 3. 



Table 2. The density and mass thresholds used to produce cluster 
catalogues with mean separations d c i = 10 — 50h~ 1 Mpc. The 
mass thresholds are given for the FOF clusters determined by 



b = 0.2. 
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n c i ~ 2.5 x 10~ 5 /i 3 Mpc -3 (dei ~ 34/r 1 Mpc), respectively 
(Dalton et al. 1994, Einasto et al. 1997). 

Fig. 4 demonstrates the rms peculiar velocity of clus- 



ters with a mean separation d c i = 10 — 50/i _1 Mpc. In this 
range, the rms velocity of DENSMAX clusters is almost in- 
dependent of the density of clusters. For the clusters with a 
mean separation d c i = 30/i _1 Mpc, u rms = 510 kms -1 and 
Vrms = 580 kins -1 in the ACDM1 and rCDM models, re- 
spectively. These values are similar to the linear rms velocity 
of peaks at the radius R ~ lfe _1 Mpc. 

The rms velocity of FOF clusters decreases with cluster 
richness. For rich clusters, the rms velocity of FOF clusters 
is smaller than the rms velocity of clusters determined with 
the DENSMAX method. In the ACDM model for b = 0.2, 
the rms velocities are 500 kms -1 and 430 kms -1 , if d c i = 
10/iT 1 Mpc and d c i = 50/i _1 Mpc, respectively. For b = 
0.15 and b — 0.3, this effect is similar. For the clusters with 
dei = 30ft _1 and b = 0.2, v rms = 475 kms -1 and v rm3 — 
565 kms~ in the ACDM1 and rCDM models, respectively. 

In Table 3 we compare the rms velocity of FOF clusters, 
Vrms, with ffp(R) for the radius R = l/i _1 Mpc. We analyzed 
the rms velocity of clusters for different values of the cluster 
separation. The results are given for the clusters determined 
by b = 0.2. The rms peculiar velocities v rm3 for the clusters 
determined by b = 0.15 and b = 0.3 are similar. In the 
ACDM1 model, the rms peculiar velocity of small clusters 
is close to the linear theory expectations, while the rms pe- 
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Table 3. Comparison of the rms velocity of FOF clusters, v rms , 
with the linear theory predictions for peculiar velocities of peaks, 
Up, for the radius R = lh~ x Mpc. The results are given for the 
FOF clusters determined by b = 0.2. 



d c i 


Vrms / °~p 


Vrms / °~p 


Vrms / °~p 


(A- 1 Mpc) 


ACDMl 


tCDM 


ACDM2 


10 


0.98 


1.03 


0.99 


30 


0.93 


1.01 


0.92 


50 


0.84 


0.94 


0.88 



culiar velocity of rich clusters is smaller (~15% per cent for 
clusters with a mean intercluster separation d c i = 50/i -1 
Mpc). 

Our results are in good agreement with the results ob- 
tained by Colberg et al. (2000). They also studied the rms 
velocity of clusters in the Virgo ACDM and rCDM mod- 
els, but used a slightly different method to select clusters. 
High-density regions were located using a FOF method with 
b — 0.05 and their barycentres were considered as candidate 
cluster centers. Any candidate centre for which mass within 
1.5/i _1 Mpc exceeded the threshold mass M t was identified 
as a candidate cluster. The final cluster list was obtained by 
deleting the lower mass candidate in all pairs separated by 
less than 1.5/i _1 Mpc. The peculiar velocity of each cluster 
was defined to be the mean peculiar velocity of all the parti- 
cles within the 1.5/i _1 Mpc sphere. In this method, the size 
of the selected clusters is same for all clusters and in this 
sense, this method is similar to the DENSMAX method. 

Colberg et al. (2000) used the value M t = 3.5 x 
1O 14 /i _1 M0. For this value, the number of clusters was w 70 
in the ACDMl and rCDM models (d ct w 58/i _1 Mpc). They 
found that the rms cluster velocities derived from simula- 
tion, v s are 439 kms -1 and 535 kms -1 in the ACDMl and 
rCDM models, respectively (they did not include the disper- 
sion v\). If we use the DENSMAX method, we find that the 
rms velocity of clusters is almost independent of the number 
density of clusters and for d c i = 50/i _1 Mpc, the velocities 
are v 3 = 450 kms -1 and 549 kms -1 in the ACDMl and 
rCDM models, respectively. These values are very close to 
the values found by Colberg et al. (2000) , only slightly larger 
(» 2 per cent). This small difference is probably caused by 
the fact that in the DENSMAX method we use the smooth- 
ing length ~ 1ft -1 Mpc, which is smaller than 1.5ft _1 Mpc 
used by Colberg et al. (2000)]. For comparison, the rms ve- 
locities v s for the FOF clusters for b = 0.2 and d c i — 50/i _1 
Mpc, are 368 kms -1 and 468 kms -1 in the ACDM and 
rCDM models, respectively. 

Let us now consider the rms velocities of clusters for dif- 
ferent cluster radii. We studied the rms velocity of clusters 
with the effective radius R e f / larger than a given threshold 
radius. The clusters were ranked according to their effective 
radius and we selected N c i = (L/d c ;) 3 highest ranked clus- 
ters to produce cluster catalogues with mean separations 
10 — 50h~ Mpc. Table 4 shows the threshold radii used for 
different values of the mean cluster separation. 

Fig. 5 illustrates the dependence of 

Vrms OH 

cluster radii. 

We see that the effect of the cluster radius on v rma is similar 
to the effect of the cluster mass on v rms (compare Fig. 4 
and Fig. 5) . In the models studied, the rms velocity of small 
clusters is higher than the rms velocity of large clusters. For 
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Figure 5. The rms peculiar velocities of clusters for different 
values of the mean cluster separation, d c i . The clusters are ranked 
according to their effective radius, R e ff- The solid lines show the 
rms velocities in the ACDMl model and the dot-dashed lines in 
the tCDM model. The upper panel shows the velocities for the 
clusters defined by b = 0.15, the middle panel for the clusters 
defined by b = 0.2 and lower panel for the clusters defined by 
b = 0.3. 
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Table 4. The threshold radii Rt used to produce cluster cata- 
logues with mean separations d c i = 10 — 50/i -1 Mpc. The radii 
are given for the FOF clusters defined by b = 0.2. 





Rt (ACDM1) 


R t (tCDM) 


(ft" 1 Mpc) 


(ft- 1 Mpc) 


(ft- 1 Mpc) 


10 


0.22 


0.24 


20 


0.40 


0.40 


30 


0.58 


0.53 


40 


0.74 


0.64 


50 


0.88 


0.72 



example, in the ACDM1 model for b — 0.2, the rms velocities 
are 520 kms -1 and 450 kms -1 for clusters with d c i — 10ft -1 
Mpc and d c i = 50h~ Mpc, respectively. 



4 THE SHIFTED-MESH SCHEME 

When we started the study of cluster velocities we first used 
a PM code with a shifted-mesh scheme. In this scheme, the 
acceleration on the grid was calculated by using the approx- 
imation 

9x(i + i,J, k) = f(i,j,k) - ip{i + l,j,k), (13) 

where tp(i,j, k) is the gravitational potential on the grid. A 
similar scheme for one-dimensional systems was proposed 
by Melott (1986). The traditional two-point finite-difference 
approximation for the acceleration is given by 

g .(M-,*) = y(i - 1 ' J ' fc) - y(i + 1 ' J '* ) . (14) 

Gramann (1987) studied the evolution of an one-dimensional 
plane-wave perturbation by using a PM code with differen- 
tial operators (13) and (14) and found that the approxima- 
tion (13) leads to smaller deviations from the exact solution 
than the approximation (14). 

Fig. 6 shows the rms peculiar velocities of clusters in 
a PM simulation, where we used the approximation (13) 
for a three-dimensional system. The cosmological parame- 
ters in this simulation were chosen similar to the ACDM1 
and ACDM2 models. We investigated the evolution of 256 3 
particles on a 256 3 grid in the comoving box of size L = 
239.5ft -1 Mpc. For comparison we show the rms peculiar 
velocities in the ACDM2 model with a standard finite- 
difference approximation (14). The clusters are defined by 
the FOF method with b = 0.2. The rms cluster velocities in 
the ACDM2 model are similar to the ACDM1 model, which 
achieves a force resolution smaller than the mean particle 
separation X p (see Fig. 3 and Fig. 4). 

We see that the shifted-mesh scheme artificially boosts 
cluster velocities. The effect is particularly strong for rich 
clusters. The rms velocities are 560 km s _1 and 735 kms -1 , 
if d c i = 10ft -1 Mpc and d ct — 50/i -1 Mpc, respectively. 
This effect is probably caused by the fact that using the 
scheme (13), we calculate the components (g x ,g y ,gz) at dif- 
ferent locations on the grid and this can lead to artificial 
enhancement of particle acceleration in one dimension. This 
effect does not arise in an one-dimensional system, where 
9v = 9z = 0. 

In the ACDM2 model we used the standard two-point 



800 




400 - 



300 1 1 1 1 1 1 1 

20 40 60 

d d (h" 1 Mpc) 

Figure 6. The rms peculiar velocities of clusters in the PM simu- 
lation, where we used the shifted-mesh scheme (13) (dashed line) 
and in the ACDM2 model using a standard approximation (14) 
(solid line) 



finite-difference approximation (14). In this model, the rms 
velocities of clusters are similar to the ACDM1 model. 
Therefore, we can use the PM code to study peculiar veloc- 
ities on scales that are close to the mean particle separation 
and force resolution. But it is important to use a correct 
difference operator in simulations. 



5 SUMMARY AND DISCUSSION 

In this paper we have examined the rms peculiar velocities of 
galaxy clusters, v rms , for different cluster masses and radii. 
We analyzed clusters in the Virgo simulations for two cos- 
mological models, ACDM and rCDM (Jenkins et al. 1998). 
These simulations were carried out using the AP 3 M code. 
We used the simulations where the mean particle separation 
\ p ~ Ih^ 1 Mpc. We also analyzed clusters in an N-body sim- 
ulation where the evolution was followed using a PM code 
with the same mass resolution that was used in the Virgo 
simulations. The cosmological parameters for this simulation 
were chosen similar to the Virgo ACDM model. We found 
that the rms velocities of clusters in the PM simulation are 
similar to the rms velocities in the AP 3 M simulation. We 
can use the PM code to study peculiar velocities on scales 
that are close to the mean particle separation X p . 

To identify clusters in the simulations we used two 
methods: the standard friends-of-friends (FOF) method and 
the method where the clusters are defined as maxima of the 
density field smoothed on the scale R ~ lftT 1 Mpc (DENS- 
MAX). The velocity of DENSMAX clusters was determined 
using the same smoothing scale as for the density field. If 
we use the DENSMAX method, the size of the selected clus- 
ters is similar for all clusters. The velocity of FOF clusters 
was defined to be the mean velocity of all the particles in the 
cluster. If we use the FOF method, the size of the high-mass 
clusters is larger than the size of the low-mass clusters. We 
studied the velocities of FOF clusters defined by the param- 
eters b — 0.15, b = 0.2 and b = 0.3 (the parameter b defines 
the linking length as bX p ). 
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Figure 7. The cluster peculiar velocity defined using the DENS- 
MAX method vs the cluster peculiar velocity defined by the 
FOF method. We show the velocities for rich clusters for d c i = 
50/i — x Mpc. The upper panel shows the velocities in the ACDM 
model and the lower panel in the rCDM model. 



we found that the fraction of matched clusters is F = 0.72. 
For d c i = 30/i _1 Mpc in both methods, the fraction F = 0.74. 
In the rCDM model, the fraction F was similar. We se- 
lect the same objects by using different methods, but we 
rank them in a somewhat different way. In the DENSMAX 
method, the clusters are ranked according to their density 
and in the FOF method according to their mass. If we com- 
pared DENSMAX clusters for d c i = 50/i _1 Mpc with the 
FOF clusters for d cl = 40/i _1 Mpc, we found that F = 0.98. 

But we do assign different velocities for the same clus- 
ters by using different methods. Fig. 7 shows the peculiar 
velocities determined with different methods for the same 
clusters. We show the velocities of rich clusters that match 
for d c i = 50/i _1 Mpc in both methods. We see that the 
velocities of clusters defined by the FOF method are sys- 
tematically smaller than the velocities of clusters defined by 
the DENSMAX method. This is probably caused by the fact 
that in the FOF method, the size of rich clusters is larger 
than the smoothing scale R ~ lfr -1 Mpc that was used to 
define the cluster velocities in the DENSMAX method. 

We compared the rms velocities of clusters with the 
linear theory predictions for the rms peculiar velocities of 
peaks, <j p , for the smoothing radius R = l/i -1 Mpc. At 
this radius, a p = 509 kms -1 and a p = 562 km s _1 in 
the ACDM and rCDM models, respectively. We analyzed 
the rms velocity of FOF clusters for different values of the 
cluster separation. In the ACDM model, the rms peculiar 
velocity of small clusters is close to the linear theory ex- 
pectations, while the rms peculiar velocity of rich clusters 
is smaller («15 per cent for clusters with a mean separa- 
tion d cl = 50/i _1 Mpc). The rms velocity of DENSMAX 
clusters is almost independent of the cluster density and 
is similar to the linear theory expectations. For the DENS- 
MAX clusters with a mean separation d c i — SOh^ 1 Mpc, 
Vrms = 510 kms -1 and v rrn s = 580 kms -1 in the ACDM 
and tCDM models, respectively. On scales probed by galaxy 
clusters, velocity fluctuations are in the quasi-linear regime. 



We found that the rms velocity of clusters defined with 
the DENSMAX method is almost independent of the density 
of clusters. The rms velocity of FOF clusters decreases with 
the cluster mass and radius. The effect of the cluster radius 
on Vrms is similar to the effect of the cluster mass on v rms 



(see Fij;. 4 and Fig. 5). The rms velocity of small clusters 



is higher than tne rms velocity ot large massive clusters. 
For different values of b, this effect is similar. In the ACDM 
model, the rms peculiar velocity of massive clusters with an 
intercluster separation d c i = 50/i -1 Mpc is » 15% smaller 
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than tl e rms velocity of clusters with a separation d c i — 

10/T 1 Mpc. 

The rms velocity of massive FOF clusters is smaller 
than the rms velocity of high-density clusters determined 
with the DENSMAX method. What is the reason for this 
difference? Do we select different objects by using different 
methods or do we define different velocities for the same 
clusters? We analyzed the fraction of DENSMAX clusters, 
F, which match FOF clusters (in terms of their positions). 
We studied the FOF clusters determined with b — 0.2. In 
the ACDM model, if we compared DENSMAX clusters for 
d c i = 50ft" 1 Mpc with the FOF clusters for d d = 50/i _1 Mpc, 



programs used in this pa per are available at http://www- 



cc.astro.washington.edu. These programs were developed 



in the University of Washington. 
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